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1.0 Intoduction and Summary 


A numbar o-f altarnativa approach** to r*nd*zvous and proximity opara- 
tion* activitias hava baan proposad, and in soma casas damonstratad, 
in various mannad spaca-f light programs from Prcjact Gamini to tha 
Spaca Shuttla. Thasa oparations hava ganarally baan limitad in scopa, 
and fraquantly dapand haavily on ground support facilitias and parson- 
nal for both planning and axacution. Tha difficulty of automating thasa 
activitias, and tha rasulting systam parformanca raquiramants, ara 
strongly influancad by tha oparational tachniquas amployad, and by tha 
ranga of prospactiva oparations to ba carriad out. 

An axamination of tha Orbital Manauvaring Vahicla Phasa B dasign rafar- 
enca missions ravaals that as a group, thay hava tha folloaing charac- 
tsristics and implications for automatad systams and oparationsi 

a. Thay ancompass a broad spactrua of oparational activitias 
ranging from sarvicing and rapair of Iom aarth orbital plat- 
forms, to ratriaval and amplacamant of payloads at altitudas 
of up to 1000 nautical milasi 

b. Nina out of ten missions involve at laast on* rendezvous, 
and seven out of ten involve two or more; 

c. Mission design techniques must take explicit account of 
differential nodal regression effects betwaan the point nf 
origination, target, and point of termination; 

d. Mission operational techniques must be capable of hand- 
day-of-mission dispersions away from nominal conditions 
without radical revision of the operational plan; 

e. A significant possibility exists that large software de- 
velopment costs will be incurred for design of phase-speci- 
fic guidance and control functions, unless a unified ap- 
proach to guidance and control can ba devised; 

As a result of these considerations, the contractor has developed a 
generalized mission design scheme which utilizes a standard mission 
profile for all OMV rendezvous oparations, recognizes typical opara- 
tional constraints, and minimizes propellant penalties due to nodal 
regression effects. This scheme has been used to demonstrata a unified 
guidance and navigation maneuver processor (tha UMP) , which supports 
all mission phases through station-keeping. Section 2.0 describes tha 
approach to OMV oparations analysis; a description of the rasulting 
mission design systam is contained in Section 3.0, together with an 
example of its application to the Large Observatory Servicing Mission, 
which is typical of OMV rendezvous oparations. The initial demonstra- 
tion version of tha Orbital Rendezvous Mission Planner (ORMP) has baan 
provided to tha customer for evaluation purposes, and program operation 
will also be briefly discussed. 

Changes and additions to tha BETEL6EU8E basic simulation program hava 
incorporated the Unified Maneuver Processor (UMP) , tha derivation and 
implementation of which ia contained in Section 4 and Appendix A. Both 
the ORMP and UMP emphasize tha pertinence of the opening statement of 
this section that the choice of oparational technique is critical to 
tha automation problem, especially as regards on-board software over- 
head and associated development costs. 
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2.0 Rvquiramsnts Analymis 

Tabl* 2-1, takan -from tha OMV Tachnlcal Raquiramants Documant (1), aua- 
aarizaa tha aarly OMV daaign ra-faranca fliisBiona. As statad aarliar, of 

thasa tan mlBsions, all 
but ona involva at laast 
a final randazvcxis Nith 
tha Orbitar, and savan 
raquira an initial ran- 
dazvouB Mith a platform. 
In addition to tha axpli- 
cit plana changa raquira- 
mants, vary larga out-of- 
plana componants arisa at 
soma point in missions 1, 
2, and 6 dua to diffaran- 
tial nodal ragrassion. As 
much as 1000 or mora faat 
par sacond is axhibitad 
in Mission 6. It is also 
notad that thasa missions 
all raturn to thair ini- 
tial starting point, tha 
Orbitar, which implias 
that tha sacond randaz- 
vous dasign is constrain- 
ed by tha targat orbit, 
and tha daploymant orbit 
of tha shuttla. It can ba 
anticipatad that unlass 
a systam such as GPS is 
avai labia, larga naviga- 
tion arrors will occur batwaan major burns, and that avan with GPS, ap- 
plication errors will causa significant dispersions at tha targeting 
end points. In either case, it seams that as many as three mission 
phases may be required: 

i> Initial boost, gross phasing, and out-of-plana correction 
to an offset point in phase and altitude! 

ii) Final boost, phasing, and out-of-plana corrections to ra- 
ramova tha major affects of accumulated navigation errors and 
dispersions! 

iii) Terminal area operations aided by vehicle— to vehicle na- 
vigation, leading to the conclusion of the rendezvous phase. 

These mission phases are fairly generic to all rendezvous operations, 
dictated by the physical limitations of on-board actuators and sansors. 
They do not of themsalvas prescribe a particular approach to mission 
dasign and operations. Additional raquiramants which must be taken into 
account in developing an approach to automated operations fall into 
three general categories: 1) those dictated by the automation objective 
itself! 2) those enforced by a desire for consumables economy and mini- 
mal demands on sensors and subsystems! and 3) safety and operability 
issues in tha terminal proximity operations zone. Tha implications of 


MISSION TITLE 

RENDEZVOUS 
PAYLOAD 1 ORB I TER 

Larga Observatory 
(1) 

350 X 350 i 160 X 160 

Payload Placement 
(2) 

: 160 X 160 

Payload Retrieval 
(3) 

37G X 378 S 160 x 160 
1 deo o/c. WTR launch 

Payload Reboost 
(4) 

260 X 260 \ 160 x 160 
inolane reboost to 370 

Payload Deboost 
(S) 

160 X 160 ! 160 X 160 
deboost to entry 

Payload Viewing 
(6) 

1400 X 1400 S 160 X 160 
1 deo o/c each wav 

Subsatellite 

<7) 

: 160 X 160 

Multiple Payload 
(B) 

270 X 270 : 160 x 160 
deolov oavload at 245 

Upper Stage 
(9) 

□MV derivative stage 
no retrieval reouired 

Early Limited 
Servicino <10) 

as required! 160 x 160 
within capabilities 


Table 2-1: Initial OMV Dasign Reference 
Missions 


(1) RFP B-1-4-PP-02142, 3 January 19B4, Attachment A, pp. 21-23 



thM« constraint catsgoriss Mill ba sunMnarizsd in rsvsraa ordsr, as 
thay tand to dat ar mi na tha functional raquiramants on tha automatad 
systam, which is tha'ultimata topic of this saction. 

SAFETY AND OPERABILITY 

-Control of tarminal phasa approach diraction and braking off- 
sat, so as to protact payloads from accidantal collision and 
pluma impingamant, and facilitata tha acquisition of station- 
kaaping and docking sansors. 

-Control of tarminal phasa arrival tima so as to parmit visual 
monitoring of tha final approach and subsaquant proximity op- 
arations activity. 

-Control of tha manauvar timas so as to assura tha oparability 
of optical sansors for navigation support, and permit moni- 
toring of burn computations and applications. 

-Graceful performance degradation in tha face of off-nominal 
initial conditions and sensor failures. 

ECONOMY AND PERFORMANCE REQUIREMENTS 

-Placement of maneuvers at or near tha optimal timas for maxi- 
mum efficiency, particularly at nodal crossing points. 

-Predictability of the pre-terminal phase conditions so as to 
limit the tracking range requirements of RF sensors. 

AUTOMATION REQUIREMENTS 

-Capability to handle off-nominal dispersions without redesign 

of the maneuver plan. 

-Feasible software implementation requirements. 

It is not necessary to search any farther than the Apollo/Sky lab con- 
centric flight plan technique to find an operational stratagem which 
meets all of the above requirements, with the exception of automated 
stationkaaping. This will ba dealt with later, after an exposition of 
the design technique afforded by CFP, and the resulting approach to au- 
mated mission design and automated orbital operations. Since it is the 
terminal conditions of the operation that are ultimately of interest, 
we begin with tha design of the terminal phase, and proceed backward to 
the initial conditions. 

TERMINAL PHASE 

Under a previous contract (2), the conditions for a minimum propellant 
transfer to a coalliptic offset from a target have bean derived. For 
the special case where tha initial orbit is coalliptic with tha target 
orbit, it can be shown that the sum of the thrust directions at trans- 
fer and braking add up to ISO degrees, in the target local vertical 
system. Furthermore, if tha final approach to the target is from an 
elevation of -26.72 degrees as seen from the target, the transfer ma- 
neuver Mill ba along the current 1 ine-of-sight at TPI, and the trans- 


(2) Multi-body Proximity Operations GNtiC Analysis Program, Software Re- 
quirements Document, Appendix B, Contract NA89-16896, 1 Dec 1983. 
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fer interval is 151.75 degrees orbit travel. These conditions are 
independent o-f orbital altitude or initial range, providing the trans- 
-fer is accomplished upon the appearance o-f the 26.72 degree elevation. 
Since it is desireabl*; to minimize excursions a.^ay from the track at- 
titude -for per-forming t^-T-ficient maneuvers, it seems reasonable to a- 
dopt these as the nominal conditions for a terminal phase. If the time 
of arrival is fixed to be shortly after local sunrise, and the rendez- 
vous is from below, the target will be illuminated from behind the in- 
'lerceptor during the final braking phase. An entirely serendipitous 
benefit of this choice results from the fact that in low earth orbit, 
the transfer interval is almost exactly the maximum duration of a 
night side pass, so that target sunset occurs at approximately the time 
of t ansfer, which is the optimum sun illumination angle for optical 
tracking in reflected sunlight. 

At any given target orbit altitude, the following events can now be 
fixed in the prospective timeline: 

1) TTSR, the time of target sunrise 

2) TTPF, the time of braking (=TTSR + DTSR) 

3) TTPI, the time of transfer, 151.75 deg of target orbit travel 
before TPF 

The targeting can aim the vehicle to any desired offset in altitude and 
downrange to avoid plume i '^pi ngement , and should include several nomi- 
nally zero midcourse corrections between TPI and TPF. 


PRE-TP I TARGETING 

In general, it is desireable to guarantee a minimum length of time 
prior to TPI within which both optical and RF (or other) tracking can 
be accomplished. The first requirement is taken care of by the previous 
specification of the TPI conditions, which guarantee target visibility 
with back illumination immediately prior to TPI. In order to accomodate 
the second requirement, we specify the minimum interval prior to TPI 
during which the range to the target must be less than the maximum 
tracking range of the RF sensor. Since the transfer optimality result 
used in specifying the transfer phase conditions assumed transfer frcm 
a coelliptical orbit, we require that the minimum dwell in the pre-TPI 
coelliptical phase be given as DTCQE, and that the maximum range at the 
onset of this interval be given as RMAX. We translate the dwell time 
into an equivalent arc of orbit travel , WTCOE, and find that the coel- 
liptical phase must begin at least (3/2) -»WTCOE-»DELTAH in curvilinear 
distance from the TPI point, since the catch-up rate in coelliptical 
phase is well known as - (3/2) -»W-»DELTAH. Since the maximum range at this 
point has been given, the required values of DELTAH and downrange dis- 
distance, DELTAS, can be solved for. We assume, for the moment, that 
the maximum range point corresponds to a cuelliptical maneuver. We have 
thus determined the additional constraints 

4) DELTAH at the coelliptical burn 

5) DELTAS at the coelliptical burn 

6) TIGCOE latest time of the coelliptical burn 

7) Four additional conditions, imposed by the requirement to be 
in— plane and coelliptical 


PRE-COELLIPTICAL TARGETING 

Conditions 4-7 above completely specify the conditions at the onset of 
the coelliptical phase. Further specification of the placement and tim- 
ing of prior maneuvers depends on whether, from a planning point of 
view, the initial state of the interceptor can be considered uncon- 
strained in some degree. The most readily available degree of freedom 
at the onset of a rend'''zvous sequence is the initial phase between the 
target and the interceptor. Since typical operations of an OMV deployed 
from either the shuttle or space station fix the orientation and shape 
of the interceptor orbit, we consider that the OMV possesses sufficient 
energy relative to the interceptor to select the initial phasing by 
varying the time of the initial maneuver in the pre-coelliptical phase. 
The objective of the pre-coelliptical program must therefor be to deal 
with the existing dispersions, differential altitude, and out-of-plane 
conditions between it and the target. Operational and system perfor — 
mance constraints inevitably set a limit on the maximum time allowed 
for the rendezvous phase of the mission. We take this limit, TMAX, as 
given, and subtract from it the time required for the phases from the 
coelliptical to the braking maneuver, TTPF-TIGCOE. In the error-free 
case, the most efficient sequence of maneuvers for creating the requir — 
ed conditions at TIGCOE result from a phasing burn at the largest in- 
tegral multiple of interceptor orbits prior to TIGCOE which is greater 
than time zero, followed by a height adjust maneuver 1/2 rev before 
TIGCOE. At the same time, it must be noted that for large differential 
altitudes, incrensinq the total time for the rendezvous sequence will 
also increase the effects of differential nodal regression, which may 
either decrease or augment the out-of-plane condition at TIGCOE. The 
first candidate solution is obtained by noting that TIGCOE is ultimate- 
ly tied to the times of target sunrise, which repeat each target orbit 
period, and choosing the latest TIGCOE which does not violate the total 
time constraint: 

1) Choose latest TIGCOE sach that TTPF .LE. TMAX; 

2) Find the first time of nodal crossing ,TNODE, prior to TIGCOE; 

3) Augment the time in the coelliptical phase by DT=TIGCOE-TNODE , 
and the downrange target by (3/2) *DELTAH*DT*W; 

4) Redefine TIGCOE=TNODE; 

5) Find the largest integral number of interceptor orbits prior to 
TIGCOE which is larger than time zero, TIGl; 

6) Define TIGNH as being 1/2 rev prior to TIGCOE: 

7) Solve for the targets with burns at TIGl, TIGNH and TIGCOE. 

The next candidate solution is obtained by repeating this process using 
the next earlier TIGCOE (one target orbit revolution), until a mimimum 
of the total propellant for the three burns is found, or until no fur- 
ther reductions in total duration are possible. This procedure guaran- 
tees a minimum solution for the total problem under given specification 
of the initial interceptor orbit, target orbit, and total time allowed. 

Further reductions in total propellant may be achieved some of the ini- 
tial interceptor orbit parameters are free variables in the design pro- 
lem. In such cases, the initial inclination or ascending node may be 
chosen so as to counteract the effects of nodal regression, or the ar — 
ument of perigee and true anomaly fixed so as to place TNODE at the 
original value of TIGCOE, and place TIGl, TIGNH and TIGCOE at relative 
apsides of the interceptor and target orbits. 

Once a minimum propellant solution has been determined, a final modifi— 
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cation IS appropriate if a priori statistics are available on the 
pected envelope o-f dispersions, defined as the difference between ac- 
tual and desired initial condition on the day of execution. The covar- 
l ance of dispersions should be propagated to TIGCOE, and all maneuvers 
forward in time by an interval which decreases the initial phase by 3- 
siqma of the expected up— range dispersion at TIGCOE. This guarantees 
that any changes in maneuvers due to actual dispersions are in the 
direction of decreasing relative energy between the initial and final 
orbits. Also, the mission timeline must provide for a (nominally zero) 
out-of-plane burn at 1/4 rev before TIGCOE to handle out-of-plane dis- 
persions on the day of execution. 


3.0 Orbital Rendezvous Mission Planner (ORMP) 


Figure 3-1 presents the input data 
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FIGURE 3-1: Planner Input Data 


specifies that the next rendezvous 


file for the ORMP, which is identi- 
cal in format to that for the 
ORBNAV variant of Program BETEL- 
GEUSE delivered under this con- 
TRACT. This program was developed 
at SSI to demonstrate the feasi- 
bility of automated mission de- 
sign using an approach proven on 
programs previous to the Shuttle. 
Of interest in Figure 3-1 are i- 
tems in the maneuver command file 
which IS set up to construct a 
two-rendezvous sequence for OMV 
OMV Design Reference Mission 1. 
On input to the ORMP, the sign of 
the variable SGN DH is examined 
to determine whether the rendez- 
vous is to be from below (0) or 
above (>0) , and the the minimum 
range at TIGCOE is given as R-NSR 
= 12 nautical miles. On cards 4 
and 5, the difference between 
TTIG and TCOE (=1800) is taken to 
be the minimL.n time to be spent 
in the coelliptic phase prior to 
TPI. H-OFF and S-OFF on card 5 
specify that the TPI targeting 
will aim for .002 n.m. uptrack of 
the target and along v-har. The 

variable DWT on cards 6, 7, and 8 
specify degrees of orbit travel 
at which each of these maneuvers 
occurs following the the one pre- 
vious. The braking maneuver of 

this particular example therefore 
occurs 139 degrees of orbit tra- 
vel after TPI. H-OFF and S-OFF 

specify the same targeting for 
the midcourses as for TPI. On 
On card 8, DTSR specifics that 
trakinq will occur .05 hours af- 
ttr target sunrise, and DWELL 

cannoi. be initiated earlier than 1 
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Figure Z~Zr Planner Output Data 


hcur after braking o-f the first 
seguence. The same definitions 
apply to the same variables of 
cards 9-16 for the second rendez- 
vous, which IS a return of the 
observatory to the initial orbit 
of the OliV. TTIG of card 16 spe- 
cifies that the total sequence 
must be CEjncluded by 27 hours af- 
ter burn 1 of the first rendez- 
vous. No specification of this 
limit on card 8 of the first ren- 
dezvous IS interpreted by the 
□RMP as freeing the argument of 
perigee and true anomaly of the 
□liV to be chosen so as to min*. - 
mize propel l:ii t reguirements as 
described above. Also included in 
Fiqt'fe 3-1 are the input initial 
initial states, covariance of 
target relative OMV dispersions, 
and covariance of OMV-relative 
target orbit uncertai nt 1 es. Fig- 
ure 3-2 exhibits the output file 
where the input values of time, 
coelliptical targets, and initial 
orbits have been replaced by 
those resulting from the ORMP 
construction of the mission. This 
This data file may be input to 
□RBNAV without change for subse- 
quent analysis of the total pro- 
pellant requi rements . Finally, 
Figures 3-3a and 3-3b summarize 
the nominal maneuyer solutions 
for the TIGl through TIGCOE se- 
quence. The method of computation 
of these solutions will be the 
subject of the next section. 


M* ii Ti^ .1 mi»i> .M mm* -i.n min* .m m iw i.o 
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MB ii Tit* MT4.2 mm* .M mm* m i«iiii* .w m i** .m 

MB 4i Til* S32t.! OVIIII* -.07 IVTIII* Ut.H Willi* -Ilf.W i« W* M2. 4] 

M* 2.42 NF* -10.71 ItF* .M WOI* .11 WCW* .M WCB* .M 

! craa MV* 41. MW* .M Ml » M.n 

Figure 3-3a: First Rendezvous 
Phasing, NCC , OOP, Coelliptic 
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Figure 3-3b; Second Rendezvous 
Phasing, NCC, OOP, Coelliptic 
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A.O Universal Manecver Processor (UMP) 

The classical solution to the minimum 4uel n-burn rendezvous problem 
has been treated by Lawden (3) and others, and shown to be a two-point 
boundary value problem. In Lawden's treatment , a primer vertor is de- 
fined which can be shown to follow the Cl ohessv-W? 1 tshi r e bquations of 
relative motion. For an impulsive sequence, it is known that the opti- 
mal maneuver points occur at time when the aapnitude of the primer vec- 
tor approaches 1 from below (<1), and that tne thrust direction is that 
of the primer vector at each such point. The boundary value problem is 
that of finding a solution to the primer vector equation such that the 
sequence of impulses result in attainment of the specified end conoi- 
tions, given specified initial conditions. In general this problem can- 
not be solved in closed form, and a different approach is more practi- 
cal . By adopting coelliptical conditions as the nominal final state of 
various phases of the rendezvous sequence, it can be calculated what 
initial phasing results in mimimal total propellant for a two-impulse 
transfer to a coelliptical final condition from a<^bitrary initial con- 
ditions. If the relative state at each point is given in terms of posi- 
tion, and velocity relative to coelliptical speed at that point, the 
resulting problem can be minimized with resnect to transfer interval, 
which amounts to a solution to the primer vector boundary value prob- 
lem. This has been done in the ORMP. 


An important feature rj-f the resulting mission plan is that dispersions 
in relative motion due to initial placement errors and resolved naviga- 
tion uncertainties dr- not significantly affect the time placement of 
maneuvers for the optimal sequence. Such effects as occu*^ can be dealt 
with by minor adjustments of the maneuver times, which do not alter the 
basic timeline. Therefore, it is almost never necessary to re-sol ve the 
planning problem in real time. Rather , it is only necessary to solve 
for a specific set of maneuvers based on the current best knowledge of 
the state. 


From an analytical point of view then, the maneuver computation problem 
on a given day of execution can be character i zed by a fixed number and 
time of burns, with specified initial and terminal conditions. This 
IS true of both the terminal phase and pr e-coel 1 i pt i c phases, and since 
in fact the terminal phase is a special case of the pre-coel 1 i pt i c 
phase resulting in coelliptic conditions at braking, all the optimiza- 
tion results derived for transfer to a coelliptic offset apply equally 
to the terminal phase. It is therefore possible to devise a single com- 
putational algorithm which accepts specification of the burn times and 
terminal constraints, and produces a minimum variance estimate of the 
burn solution set. Such an algorithm has been devised as described in 
Appendix A of this report, and incorporated into both the ORNP and the 
□RPNAV programs. The algorithm possesses the following important pro- 
per t i es: 


a) No computat 1 onal singularities, provided that a solution for 
a given set of times and terminal constraints physically exists; 

b) Spec 1 f i c at 1 (jn of 1 to 6 terminal constraints in position and 
velocity; 


(3) Lawden, D.F., Optimal Trajectories for Space Navigation, Butter- 
worths, London, 1963 
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c) Speci -f 1 cat i on of velocity offsets with respect to in-plarie 
coelliptical speed; 

d) Specification of out-of-plane position and velocity off-sets; 

e) Minimizes the squared magnitude of each component of each 
maneuver, and hence the sum squares of the total sequence; 

f) Utilizes software functions typically provided fo?* state and 
covariance propagation, and filter covariance measurement up- 
dating; 

q) Implicitly accounts for all modelled disturbing accelerations. 

It should be emphasized that these character i sti cs do not depend on any 
assumption of coelliptical initial or final conditions. The coellipti- 
cal conditions are utilized only for mission design purposes to find 
maneuver times and targeting conditions which are at or near the mini- 
mum total propellant points that would be identified by an exact solu- 
tion of the Lawden primer vector boundary value problem. The UMP can be 
used to find the minimum-propellant solution to any set of fixed burn 
times which have been derived by any other method, whether or not known 
to be optimal, such as the current shuttle rendezvous baseline. 

It was stated in Section 2 that the problem of station-keeping remained 
to be dealt with. The application of the UMP in this context is immedi- 
ate and obvious; One has only to specify the desired stationkeeping 
position, and periodically cycle the UMP on to compute a two-impulse 
solution from the current tn the desired position. Other approaches are 
possible but depend on specific and as yet undefined applications. 
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GENERAL SOLUTION TO TOE FIXED-TIME MANEUVER 
COMPUTATION PROBLEM 


1. Mathematical Preliminaries 

We are concerned with finding the best estimate of various vector-valued random 
variables. In order to facilitate the analysis, the method of ortogonal projec- 
tion will be employed. Using the Dirac notation, we define 


1 X> * X 

a vector-valued random variable 

<y| “ 1 

the adjoint of |y> 

T 

<y|x> - X = 


|xXy| - E( 2 ^^) 

T 

the expectation of ^ 


Further, the orthogonal complement of jx> with respect to jy> is defined as 

I x> * |x> - A|y> 

with A chosen so as to mimimize J * |xXxl ; 

J = [|x> - Ajy>][<x| - <y|A^] 

in the usual way, we take the variation in J with respect to a variatir-i A 
and require the result to be identically zero for an extremum* 

6J =-[|x> - A|y>]<yjo'J - 6A|y>[<x| - <y|A^] 

s (If 

which implies 

[|x> - Ajy>]<y| = |\Xyi - A|yXy| 

-1 

A * I xXy| lyXy 1 

so that 

|x> - |x> - |xXy| lyXyl" |y> 

- |x> - |x> 

|x> is called the orthogonal projection of |x> onto |y>, or in more mundane 
terms, the best estimate of |x> given jy>. It is straigliL forward to establish 
the following further results: 

|xXx| 

|xXy| 
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|xXx| = jxXxl “ jxXx] 

|xXx| = |xXx| = |xXx| 

As a practical example of the application of these results > consider an ob- 
servation process where 

|y> = H[x> + |n>; |xXx| = P, |nXn| « R, [xXn] 

If there is no a priori estimate of |x>, the best estimate of |x> is given by 
|x> = jxXy] |yXy|"^ly> 

= ph^[hph^ + 

a familiar result. If an estimate of |x> has been obtained previously, we 
construct a new estimate as a linear combination of the old estimate and the 
observation |y>. This is done by subtracting the projection of |x> onto |y> 
such that 

|y> = |y> - H|x> 
which has the property that 

lyX^I = [|y> - Hrx>]<^| 

= [H|x> - H|x>]<x| 

= H|xXx| 

= <> 

since by definition |x> = |x> - |x> and it has been shown that |xXx| =0 . 
Then 

|x'^> = |x> + lxXy|lyXy| ^|y> 

= |x> + PH^[HPH^ + R]’^|y> 

Also, by direct calculation 

P"^ = lx^><x‘^| = P - PH^[HPH^ + R]'^HP 
2.0 Fixed Time Maneuve r Sequep-.es 

Let |>f> be the relative state of the interceptor at tf , and be the 
initial state. Define 
- 

'^fo ®fo 
Afo Bfo 
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Bfn 

®fn 


Then for any sequence of n maneuvers |v^> at times t^, the final state is 
given by 






Further define 


d> = X.> - X > 

f fo o 


F - G^2 


1V>= 


vi^' 

y2> 


Then 

and further, 


- - ^fn^ 


d> = F|V> 


|V> = IVXdl IdXd 


-1 


d> 


If no a priori statistics are available on the space of V>, we are at li- 
berty to select the gauge condition |VXV| = I, so that dXdl=F|V><V| F^=FF^ 


|V> = F'^[tT’^]“V> 


which is the unweighted least squares estimate of |V>. Let U=|VXV|^^ith an 
initial value U=I. Then we may at each point define |d> to be the i ^ com- 
ponent of the deviation vector l^f^ ■ ^fol^o^’ where |X^> is the desired 

final state of the system. Thu' it is possible to recursively process the 
constraint equation: 


= |v^> + u^h^[h!Ujh.]-i|di> 
‘'i.H = >'i - UitlJhTu^b. ]-lh>. 


i= 1, m; m £ 6 
|Vj> = |0> 



If the we define D = |dXd|, then the covariance of the solution vector (not 
the error in the solution vector) is given by 

D = lvxv| = f'^Lff^]'^d[ff'^]‘^f 
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Ivxvl 




2.1 Conditions for Existence of 



It is apparent that |V> exists if the inverse [FF^]'^ 


exists, this will be 


T 

true if the matrix FT has full rank, i.e. its determinant is not zero. We 
need only examine the deterministic case, that is the case where the number 
of maneuver components is equal to the number of terminal constraints. If 
the number of components is less than the num' 2 r of constraint equations, a 
solution will not in general exist except for unique values of |d>. If we 
assume two three-asxis bums at different times, the F matrix looks like 



1 

CO 

2(l-CfP 

0 

Sf2 

2(1-c^2) 

0 


-2(l-c^^) 


0 

-2(1-Cf2) 

4s^2"2wt£2 

0 

1 

0 

0 

^fl 

0 

0 

Sf2 

w 

w(2-Cfj^) 

w(3wt^^-2s^^ 

) 0 

w(2-c^2) 

w(3wt^2“^s^2) 

0 


-wsfi 

w(2Cfi-l) 

0 

-wSf2 

w(2Cf2-l) 

0 


0 

0 

-WCfi 

0 

0 

-wcf2 


Lxamining the out-of -plane submatrix, we see that its determinant is 


w(Sf2Cfi 


SfiCf 2 ) = w sin(wt 2 ^) 


This will be zero for n(pi), n=0,l,2... It is further apparent that if 

wt^^= wt ^2 2n(pi), that column 1 of F has the same value as coulmn 4, and 

that F is therefore singular. We therefore find out-of-plane singularities if 
all bums are n(pi) apart, and in-plane singularities if all bums are 2n(pi) 
apart. This condition can be observed in the recursive solution process by 
testing ry 

h!u.h. > 0 

-1 1-1 


whereby if it fails, processing of that component should be aborted, and an 
indicator set to draw attention to the improper specification of the bum 
intervals. 
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